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Abstract 

In our work on the development of a model-independent data analysis method for esti- 
mating the spin-independent (SI) scalar coupling of Weakly Interacting Massive Particles 
(WIMPs) on nucleons by using measured recoil energies from direct Dark Matter detec- 
tion experiments directly, it was assumed that the analyzed data sets are background-free, 
i.e., all events are WIMP signals. In this article, as a more realistic study, we take into 
account a fraction of possible residue background events, which pass all discrimination cri- 
teria and then mix with other real WIMP-induced events in our data sets. Our simulations 
show that, for the estimation of the SI WIMP-nucleon coupling, the maximal acceptable 
fraction of residue background events in the analyzed data set of 0(50) total events is 
~ 10% - 20%. For a WIMP mass of 100 GeV and 20% residue background events, the 
systematic deviation of the reconstructed SI WIMP coupling (with a reconstructed WIMP 
mass) would in principle be ~ -|-13% with a statistical uncertainty of ~1^4% (~ ~3.3%^|3y° 
for background-free data sets). 



1 Introduction 



Currently, direct Dark Matter detection experiments searching for Weakly Interacting Massive 
Particles (WIMPs) are one of the promising methods for understanding the nature of Dark 
Matter (DM) and identifying them among new particles produced at colliders as well as recon- 
structing the (sub) structure of our Galactic halo [1, 2, 3, 4]. To this aim, model-independent 
methods for determining the WIMP mass [5, 6] as well as estimating the spin-independent (SI) 
WIMP coupling on nucleons [7, 8] from direct detection experiments have been developed^. 

These methods built basically on the work on the reconstruction of the (moments of the) one- 
dimensional velocity distribution function of halo WIMPs by using experimental data (measured 
recoil energies) directly [12]. The spectrum of recoil energy is proportional to an integral over the 
one-dimensional WIMP velocity distribution, fi{v), where v is the absolute value of the WIMP 
velocity in the laboratory frame. In fact, this integral is just the minus-first moment of the 
velocity distribution function, which can be estimated from experimental data directly [12, 6]. 
Then, by assuming that the SI WIMP-nucleus interaction dominates and the WIMP couphngs 
on protons and on neutrons are approximately equal, this SI WIMP coupling on nucleons can be 
estimated from experimental data directly [7, 8]. It was found that, by combining experimental 
data sets with different target nuclei, the SI WIMP-nucleon coupling can be estimated without 
making any assumption about the velocity distribution of halo WIMPs nor prior knowledge 
about the WIMP mass [7, 8] 2. 

In the work on the development of these model-independent data analysis procedures for ex- 
tracting WIMP properties from direct detection experiments, it was assumed that the analyzed 
data sets are background-free, i.e., all events are WIMP signals. Active background discrimi- 
nation techniques should make this condition possible. For example, the ratio of the ionization 
to recoil energy, the so-called "ionization yield" , used in the CDMS-II experiment provides an 
event-by-event rejection of electron recoil events to be better than 10^^ misidentification [18]. 
By combining the "phonon pulse timing parameter", the rejection ability of the misidentified 
electron recoils (most of them are "surface events" with sufficiently reduced ionization energies) 
can be improved to be < 10^^ [18]. Moreover, as demonstrated by the CRESST collabora- 
tion [19], by means of inserting a scintillating foil, which causes some additional scintillation 
light for events induced by a-dccay of ^""^^Po and thus shifts the pulse shapes of these events 
faster than pulses induced by WIMP interactions in the crystal, the pulse shape discrimination 
(PSD) technique can then easily distinguish WIMP-induced nuclear recoils from those induced 
by backgrounds'^. 

However, as the most important issue in all underground experiments, the signal identification 
ability and possible residue background events which pass all discrimination criteria and then 
mix with other real WIMP-induced events in analyzed data sets should also be considered. 
Therefore, in this article, as a more realistic study, we follow our works on the effects of residue 

"^In the literature, another method based on the maximum hkehhood analysis has also been discussed [9, 10, 11]. 
However, in contrast to the model-independent procedures, this maximum likelihood analysis requires prior 
knowledge/assumptions about the velocity distribution function and the local density of halo WIMPs. The 
WIMP mass and the SI cross section on nucleons determined by this method arc also coupled. 

^Note that, as will be discussed in more details later, the WIMP mass and the local WIMP density are needed 
for this estimation. While the former can be determined from (other) direct detection experiments directly [5, 6], 
the latter has conventionally been estimated by means of the measurement of the rotation curve of the Milky 
Way with an uncertainty of a factor of ~ 2 [3, 4]. However, some new techniques have recently been developed 
for determining the local Dark Matter density with a higher precision [13, 14, 15, 16, 17]. 

^For more details about background discrimination techniques and status in currently running and projected 
direct detection experiments see e.g., Refs. [20, 21, 22]. 



background events in direct Dark Matter detection experiments [23, 24] and want to study how 
well wc could estimate the SI WIMP-micleon coupling model-independently by using "impure" 
data sets and how "dirty" these data sets could be to be still useful. 

The remainder of this article is organized as follows. In Sec. 2 I review the model-independent 
method for estimating the SI WIMP coupling on nucleons by using experimental data sets 
directly. In Sec. 3 the effects of residue background events in the analyzed data sets on the 
measured energy spectrum as well as on the reconstructed WIMP mass will briefly be discussed. 
In Sec. 4 I show numerical results of the reconstructed SI WIMP-nucleon coupling by using 
mixed data sets with different fractions of residue background events based on Monte Carlo 
simulations. I conclude in Sec. 5. Some technical details will be given in an appendix. 

2 Method for estimating the SI WIMP— nucleon coupling 

The basic expression for the differential event rate for elastic WIMP-nucleus scattering is given 
by [3]: 



fi{v) 



dv. (1) 



Here R is the direct detection event rate, i.e., the number of events per unit time and unit mass 
of detector material, Q is the energy deposited in the detector, F{Q) is the elastic nuclear form 
factor, fi(v) is the one-dimensional velocity distribution function of the WIMPs impinging on 
the detector, v is the absolute value of the WIMP velocity in the laboratory frame. The constant 
coefficient A is defined as 

(2) 

where po is the WIMP density near the Earth and (Tq is the total cross section ignoring the form 
factor suppression. The reduced mass mr,N is defined by 

"^r,N = \ , \0) 

where is the WIMP mass and mN that of the target nucleus. Finally, Vrain is the minimal 
incoming velocity of incident WIMPs that can deposit the energy Q in the detector: 



Vmm = a^jQ, (4) 
with the transformation constant 



and i'max is the maximal WIMP velocity in the Earth's reference frame, which is related to the 
escape velocity from our Galaxy at the position of the Solar system, Vcsc ^ 600 km/s. 

For spin-independent scalar WIMP interactions, the total cross section in Eq. (2) can be 
expressed as [3, 4] 

= (^) ml^ [ZU + (A - Z)U] ' . (6) 



Here mr,N is the reduced mass defined in Eq. (3), Z is the atomic number of the target nucleus, 
i.e., the number of protons, A is the atomic mass number. A — Z is then the number of neutrons. 



/(p,n) are the effective scalar couplings of WIMPs on protons p and on neutrons n, respectively. 
Here wc have to sum over the couplings on each nucleon before squaring because the wavelength 
associated with the momentum transfer is comparable to or larger than the size of the nucleus, the 
so-called "coherence effect" . In addition, for the lightest supersymmetric ncutralino, and for all 
WIMPs which interact primarily through Higgs exchange, the scalar couplings arc approximately 
the same on protons and on neutrons: /n ~ /p. Thus the "pointhke" cross section ctq^ in Eq. (6) 
can be written as 



where mr,p is the reduced mass of the WIMP mass and the proton mass rup, and 

•4' 



7r 



<pi/pr (8) 



is the SI WlMP-nucleon cross section. Here the tiny mass difference between a proton and a 
neutron has been neglected. 

It was found that, by using a time-averaged recoil spectrum, and assuming that no direc- 
tional information exists, the normalized one-dimensional velocity distribution function of halo 
WIMPs, fi{v), can be solved from Eq. (1) analytically [12] and, consequently, its generalized 
moments can be estimated by [12, 6]*^ 



/"f (Qmax) 

K)(t'(gmin),t'(gmax)) = / f,{v) dv 

•^"(Qmin) 



2Q::^;^^ r(g^in)/F^(Qmin) + {n+ l)/n(Qmin, Q. 

1 ^2 

^Qmin^(Qmin)/ -^^(Qmin) H" -^o(Qmin7 Qmax) 



(9) 



Here v{Q) = a^/Q, (3(min,max) are the experimental minimal and maximal cut-off energies of the 
data set, respectively. 



\ ^ / expt, Q=Qmin 

is an estimated value of the measured recoil spectrum {dR/dQ)cxpt {before normalized by an 
experimental exposure S) at Q = Qmm, and /n(Qmin, Qmax) can be estimated through the sum: 

A^tot Q(n-l)/2 

-^n(Qmin) Qmax) = / , 7-19/^ \ ) (H) 

where the sum runs over all events in the data set that satisfy Qa £ [Qmin, <5max] and A^tot is the 
number of such events. Note that, firstly, by using the second line of Eq. (9) (v"')(v(Qmm), '^(Qmax)) 
can be determined independently of the local WIMP density Pq, of the velocity distribution 
function of incident WIMPs, fi{v), as well as of the WlMP-nucleus cross section (Tq. Secondly, 
f^iQmin) and In{Qmin,Qmax) are two key quantities for our analysis, which can be estimated ei- 
ther from a functional form of the recoil spectrum or from experimental data (i.e., the measured 
recoil energies) directly^. 



^Here we have implicitly assumed that Qmax is so large that a term 2Q^^^^'^r{Qma.x)/F^{Qma.x) is negligible. 
^All formulae needed for estimating r((5min), -^^(Qmin, <3max), and their statistical errors are given in the 
appendix. 



By substituting the first expression in Eq. (7) into Eq. (1), and using the fact that the integral 
over the one-dimensional WIMP velocity distribution on the right-hand side of Eq. (1) is the 
minus-first moment of this distribution, which can be estimated by Eq. (9) with n — —1, we 
have 
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(12) 



Using the definition (5) of a, the squared SI WIMP coupling on protons (nucleons) can then be 
expressed as [7, 8] 
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(13) 



Note that the experimental exposure £ appearing in the denominator relates the actual counting 
rate (dR/dQ) e^pt to the normalized rate in Eq. (1). 

As mentioned in the introduction, the Dark Matter density at the position of the Solar system, 
Po, appearing in the denominator of the expression (13) for estimating |/p|^ has conventionally 
been estimated by means of the measurement of the rotation curve of the Milky Way. The 
currently most commonly used value for po is [3, 4] 



Po ^ 0.3 GeV/cm^ 



(14) 



with an uncertainty of a factor of ~ 2. However, some new techniques have been developed for 
determining po with a higher precision [13, 14, 15, 16, 17]. These estimates give rather larger 
values for po; e.g.. Catena and UUio gave [13] 



Po = 0.39 ± 0.03 GeV/cm^ , 
and Salucci et al. even gave [15] 

Po = 0.43 ± 0.11 ± 0.10 GeV/cm^ . 



(15) 



(16) 



Moreover, instead of a spherical symmetric density profile assumed in Refs. [13, 15], in Refs. [14, 
16, 17] the authors considered an axisymmetric density profile for a fiattened Galactic Dark 
Matter halo [25] caused by the disk structure of the luminous baryonic component. It was found 
that the local density of such a non-spherical DM halo could be enhanced by ~ 20% or larger 
[14, 16] and Pato et al. gave therefore [16] 



Po = 0.466 ± 0.033(stat) ± 0.077(syst) GeV/cm^ 



(17) 



Nevertheless, since the squared SI WIMP-nucleon coupling |/pp is inversely proportional to the 
local WIMP density, by using Eq. (13) one can at least give an upper bound on |/pp. Moreover, 
as shown in Refs. [7, 8], in spite of the very few (0(50)) events from one experiment, for a 
WIMP mass of 100 GeV, the SI WIMP-nucleon coupling |/p| can be estimated with a statistical 
uncertainty of only ~ 15%; it leads to an uncertainty on the SI WIMP-nucleon cross section of 
~ 30%, which is (much) smaller than the uncertainty on the estimate of the local Dark Matter 
density. 



3 Effects of residue background events 



In this section I first show some numerical results of the energy spectrum of WIMP recoil signals 
mixed with a few background events. Then I review the effects of residue background events in 
the analyzed data sets on the reconstruction of the WIMP mass m^. 

For generating WIMP-induced signals, we use the shifted MaxweUian velocity distribution 
[2, 3, 12]: 



l,sh 



VeVQ 



(18) 



with vq ~ 220 km/s and — 1.05 vq, which are the Sun's orbital velocity and the Earth's 
velocity in the Galactic frame^, respectively; the maximal cut-off of the velocity distribution 
function has been set as Wmax = 700 km/s. The commonly used elastic nuclear form factor for 
the SI cross section [26, 3, 4]: 



{Q) 
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(19) 



will also be used^. Meanwhile, in order to check the need of a prior knowledge about an (exact) 
form of the residue background spectrum, two forms for the background spectrum have been 
considered. The simplest choice is a constant spectrum: 



(20) 




More realistically, we use the target-dependent exponential form introduced in Ref. 
residue background spectrum: 



— I = exp 

bg,ex 



g/keV \ 

^0.6 ) 



for the 



(21) 



Here Q is the recoil energy, A is the atomic mass number of the target nucleus. The power index 
of A, 0.6, is an empirical constant, which has been chosen so that the exponential background 
spectrum is somehow similar to, but still different from the expected recoil spectrum of the 
target nucleus; otherwise, there is in practice no difference between the WIMP scattering and 
background spectra. Note that, among different possible choices, we use in our simulations 
the atomic mass number A as the simplest, unique characteristic parameter in the general 
analytic form (21) for defining the residue background spectrum for different target nuclei. 
However, it does not mean that the (superposition of the real) background spectra would depend 
simply/primarily on A or on the mass of the target nucleus, ttin. In other words, it is practically 
equivalent to use expression (21) or (di?/(iQ)bg,cx ~ q-QI^"^-^ ^ev (jjrectly for a ^^Ge target. 

Note also that, firstly, as argued in Ref. [23], two forms of background spectrum given above 
arc rather naive; however, since we consider here only a few residue background events induced 
by perhaps two or more different sources, which pass all discrimination criteria, and then mix 
with other WIMP-induced events in our data sets of 0(50) total events, exact forms of different 
background spectra are actually not very important and these two spectra, in particular, the 



^The time dependence of the Earth's velocity will be ignored in our simulations. 

^Othcr commonly used analytic forms for the one- dimensional WIMP velocity distribution as well as for the 
elastic nuclear form factor for the SI WIMP-nucleus cross section can be found in Refs. [12, 11]. 



exponential one, should practically not be unrealistic^. Secondly, as demonstrated in Refs. [7, 8] 

and reviewed in the previous section, the model-independent data analysis procedure for esti- 
mating the SI WIMP-nucleon coupling requires only measured recoil energies (induced mostly by 
WIMPs and occasionally by background sources) from direct detection experiments. Therefore, 
for applying this method to future real data, a prior knowledge about (different) background 
source(s) is not required at all. 

Moreover, for our numerical simulations presented here as well as in the next section, the 
actual numbers of signal and background events in each simulated experiment are Poisson- 
distributed around their expectation values independently; and the total event number recorded 
in one experiment is then the sum of these two numbers. Additionally, we assumed that all 
experimental systematic uncertainties as well as the uncertainty on the measurement of the 
recoil energy could be ignored. The energy resolution of most existing detectors is so good that 
its error can be neglected compared to the statistical uncertainty for the foreseeable future with 
pretty few events. 

3.1 On the measured energy spectrum 

In Figs. 1 I show measured energy spectra (solid red histograms) for a ^^Ge target with six 
different WIMP masses: 10, 25, 50, 100, 250, and 500 GeV based on Monte Carlo simulations. 
The dotted blue curves are the elastic WIMP-nucleus scattering spectra, whereas the dashed 
green curves are the exponential background spectra given in Eq. (21), which have been normal- 
ized so that the ratios of the areas under these background spectra to those under the (dotted 
blue) WIMP scattering spectra are equal to the background-signal ratio in the whole data sets 
(e.g., 20% backgrounds to 80% signals shown in Figs. 1). The experimental threshold energies 
have been assumed to be neghgible and the maximal cut-off energies are set as 100 keV. The 
background windows (the possible energy ranges in which residue background events exist) have 
been assumed to be the same as the experimental possible energy ranges. 5,000 experiments 
with 500 total events on average in each experiment have been simulated. 

Remind that the measured energy spectra shown here are averaged over the simulated ex- 
periments. Five bins with linear increased bin widths have been used for binning generated 
signal and background events. As argued in Ref. [12], for reconstructing the one-dimensional 
WIMP velocity distribution function, this unusual, particular binning has been chosen in order 
to accumulate more events in high energy ranges and thus to reduce the statistical uncertainties 
in high velocity ranges. However, as shown in Sec. 2, for the estimation of the SI WIMP-nucleon 
couphng (as well as for the determination of the WIMP mass [6]), one needs either events in 
the first energy bin or all events in the whole data set. Hence, there is in practice no difference 
between using an equal bin width for all bins or a (linear) increased bin widths. 

It can be found in Figs. 1 that the shape of the WIMP scattering spectrum depends highly on 
the WIMP mass: for light WIMPs {m^ < 50 GeV), the recoil spectra drop sharply with increasing 
recoil energies, while for heavy WIMPs (m^^ > 100 GeV), the spectra become flatter. In contrast, 
the exponential background spectra shown here depend only on the target mass and are rather 
flatter (sharper) for light (heavy) WIMP masses compared to the WIMP scattering spectra. This 
means that, once input WIMPs are light (heavy), background events would contribute relatively 
more to high (low) energy ranges, and, consequently, the measured energy spectra would mimic 
scattering spectra induced by heavier (lighter) WIMPs. Moreover, for heavy WIMP masses, 
since background events would contribute relatively more to low energy ranges, the estimated 

^Other (more realistic) forms for background spectrum (perhaps also for some specified targets/experiments) 
can be tested on the AMIDAS website [27, 28]. 



'^^Ge, Q^^ < 1 00 keV, 0^,3^, bg < 1 00 500 events (20% 
AM I DAS littpj'/pisfvO.pil.physlk.iini-liiebirigen .da/daikmatler/amidas/ 



. (dR/dO]3,g„a| 

■ (dR/dQ)n,ea,^refl 



^^Ge, Qn,ax < 100 keV, Q|„a^ bg < 100 keV, 500 events (20% exponential bg), n 
AM I DAS littpj'/pisfvO.pil.physlk.iini-liiebirigen .da/darttmaHer/amidas/ 



(dR/dQ],,g^al 

(dR/dQ)t,a^g„, 

(dR/dQ)^ea^rg 



Q [keV] 



Q [keV] 



'"Ge, Q^a^ < 1 00 keV, 0^^^ bg *= 1 00 keV, 500 events (20% exponential bg), = 50 GeV 

AMIDAS httpj'/pisrvO.pit.phi'sik.jni-ljebirgEn.de/dHfkmaller/amidas/ 



'^Ge, Qma^< 100 keV, Q^aj;,bg < 100 keV, 500 events (20% exponential bg), = 100 GeV 



AMIDAS httpj'/pisrvO.pit.physik.jni-ljebmge 



(dR/dQ)„g„3| 

(dR/dQ)i,ack9rounfl 

(dR/dO)^ga5^,ga 



'> 1.4 ■ 



(dR/dQ)„g„3| 

(dR/dQ)i,ack9,oynfl 

(dR/dQ)^ga5^,ga 



10 20 30 40 50 60 70 

Q [keV] 



90 100 



10 20 



40 50 60 70 80 

Q [keV] 



'Ge, Qniax< 100 keV, Q^g^.bg < 100 keV, 500 events (20% exponential bg), = 250 GeV 



AMIDAS httpjVpisrvO.pit.physik.j 


i-Ljebmgsn de/daAmatler'am 


das/ 












(dR/dQ),ig,gi 










(dR/dQ)i,a,,g,„,,fl 








. (dR/dQ)„e3,^rea 


\ 

\ 










\ 

\ 

\ 

\ 


\ 

N 

S 










N 

•^-.^ 



















10 20 30 40 50 60 70 80 90 100 

Q [keV] 



'Ge, Qniax< 100 keV, Q^a^.bg < 100 keV, 500 events (20% exponential bg), = 500 GeV 



4.5 





AMIDAS httpjVpisrvO.pit.physik.j 


i-Ljebmgsn.de/dariimarier/anildas/ 




4 










(dR/dO),ig,«i 

(dR/dQ)t,a^g„,,^ 


3.5 










(dWdO)n,eas^rgtl 


3 












2.5 












2 












1.5 


\ 
\ 

\ 

\ 












\ 

\ 

\ 

\ 

\ 












\ 

S 








0.5 




N 





























10 20 30 40 50 60 70 80 90 100 



Q [keV] 



Figure 1: Measured energy spectra (solid red histograms) for a ^^Ge target with six different 
WIMP masses: 10, 25, 50, 100, 250, and 500 GeV. The dotted blue curves are the elastic WIMP- 
nucleus scattering spectra, whereas the dashed green curves are the exponential background 
spectra normalized to fit to the chosen background ratio, which has been set as 20% here. The 
experimental threshold energies have been assumed to be negligible and the maximal cut-off 
energies are set as 100 keV. The background windows have been assumed to be the same as the 
experimental possible energy ranges. 5,000 experiments with 500 total events on average in each 
experiment have been simulated. See the text for further details (plots from Ref. [23]). 
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Figure 2: The reconstructed WIMP masses and the lower and upper bounds of the la statistical 
uncertainties by using mixed data sets from WIMP-induced and background events as functions 
of the input WIMP mass. ^^Si and ^^Ge have been chosen as two target nuclei. The exponential 
(upper) and constant (lower) forms given in Eqs. (21) and (20) have been used for the background 
spectrum. The background ratios shown here are no background (dashed green), 5% (dotted 
magenta), 10% (long-dotted blue), 20% (solid red), and 40% (dash-dotted cyan) background 
events in the analyzed data sets in the experimental energy ranges between and 100 keV. Each 
experiment contains 50 total events on average. Other parameters are as in Figs. 1. See the text 
for further details (plots from Ref. [23]). 



value of the measured recoil spectrum at the lowest experimental cut-off energy, r{Quiin), could 
thus be (strongly) overestimated. 

More detailed illustrations and discussions about the effects of residue background events on 
the measured energy spectrum can be found in Ref. [23]. 

3.2 On the reconstructed WIMP mass 

Figs. 2 show the reconstructed WIMP masses and the lower and upper bounds of the la statistical 
uncertainties by means of the model-independent procedure introduced in Refs. [5, 6] with mixed 
data sets from WIMP-induced and background events as functions of the input WIMP mass. 
As in Ref. [6], ^^Si and ^^Ge have been chosen as two target nuclei. The experimental threshold 
energies of two experiments have been assumed to be negligible and the maximal cut-off energies 
are set the same as 100 keV. The exponential (upper) and constant (lower) forms given in 
Eqs. (21) and (20) have been used for the background spectrum. The background windows are 
set the same as the experimental possible energy ranges for both experiments. The background 
ratios shown here are no background (dashed green), 5% (dotted magenta), 10% (long-dotted 
blue), 20% (solid red), and 40% (dash-dotted cyan) background events in the analyzed data 
sets. 2 X 5,000 experiments have been simulated. Each experiment contains 50 total events 
on average. Note that all events recorded in our data sets are treated as WIMP signals in the 
analysis, although statistically we know that a fraction of these events could be backgrounds. 

Prom the upper frame of Figs. 2 it can be seen clearly that, for hght WIMP masses (m^^ < 100 
GeV), caused by the relatively flatter background spectrum (compared to the scattering spec- 
trum induced by light WIMPs), the energy spectrum of all recorded events would mimic a 
scattering spectrum induced by WIMPs with a relatively heavier mass, and, consequently, the 
reconstructed WIMP masses as well as the statistical uncertainty intervals could be overesti- 
mated. In contrast, for heavy WIMP masses (m^ ^100 GeV), caused by the relatively sharper 
background spectrum, relatively more background events contribute to low energy ranges, and 
the energy spectrum of all recorded events would mimic a scattering spectrum induced by WIMPs 
with a relatively lighter mass. Hence, the reconstructed WIMP masses as well as the statistical 
uncertainty intervals could be underestimated. 

As a comparison, the lower frame of Figs. 2 shows that, since the constant background 
spectrum is flatter for all WIMP masses^, background events contribute always relatively more 
to high energy ranges, and the measured energy spectra would thus always mimic scattering 
spectra induced by heavier WIMPs. Therefore, the reconstructed WIMP masses as well as the 
statistical uncertainty intervals are overestimated for all input WIMP masses. 

Moreover, Figs. 2 show that the larger the fraction of background events in the analyzed 
data sets, the more strongly over-/underestimated the reconstructed WIMP masses as well as 
the statistical uncertainty intervals. Nevertheless, it can be found that, with 10% - 20% 
residue background events in the analyzed data sets of ~ 50 total events, the la statistical 
uncertainty band could in principle cover the true WIMP mass pretty well. 

More detailed illustrations and discussions about the effects of residue background events on 
the determination of the WIMP mass can be found in Ref. [23] . 

^Illustrations and detailed discussions about the effects of the constant form of the residue background spec- 
trum on the measured energy spectrum for different input WIMP masses can be found in Ref. [23]. 



4 Results of the reconstructed SI WIMP— nucleon cou- 



pling 

In this section I present simulation results of the reconstructed SI WIMP coupling on nucleons^° 
by means of the model-independent method described in Sec. 2 with mixed data sets from 
WIMP-induced and background events. The WIMP mass appearing in the expression 
(13) for estimating |/p|^ has been assumed to be known precisely from other (e.g., coUider) 
experiments with an overall uncertainty of 5% of the input (true) WIMP mass or determined 
from other two direct detection experiments^^. The SI WIMP-nucleon cross section for our 
simulations is set as cr^p = 10^^ pb, the currently most commonly used value for the local 
WIMP density, po = 0.3 GcV/cm^, needed in Eq. (13) has been used for both simulations 
and data analyses. A '''^Ge nucleus has been chosen as our detector target for reconstructing 
l/pp, whereas a ^^Si target and a second ^^Ge target have been used for determining m^. The 
experimental threshold energies of all experiments have been assumed to be negligible and the 
maximal cut-off energies arc set the same as 100 keV. The exponential background spectrum 
given in Eq. (21) has been used for generating background events in windows of the entire 
experimental possible ranges. (3 x) 5,000 experiments have been simulated. 

Fig. 3 shows the reconstructed squared SI WIMP-nucleon couplings, |/p|^, and their lower 
and upper bounds of the la statistical uncertainties by using mixed data sets as functions of 
the input WIMP mass. The WIMP mass needed in Eq. (13) has been assumed to be known 
precisely. The background ratios shown here are no background (dashed green), 10% (long- 
dotted blue), 20% (sohd red), and 40% (dash-dotted cyan) background events in the analyzed 
data sets in the experimental energy ranges between and 100 keV. Each experiment contains 
50 total events on average. Remind that all events recorded in our data sets are treated as 
WIMP signals in the analysis. 

It can be found in Fig. 3 that the larger the background ratio in the analyzed data set, the 
more strongly overestimated the reconstructed SI WIMP-nucleon coupling for all input WIMP 
masses. This can be understood from the expression (1) for the different event rate dR/dQ. 
For a given WIMP mass and a specified target nucleus, the SI WIMP-nucleus cross section is 
proportional to the total event number-to-exposure ratio^^: 




Here and throughout the subscripts "sg" and "bg" stand for WIMP signals and background 
events, respectively; Sgg is the required exposure to observe the expected "WIMP signal" {not 
total) events. Note that, since background events in our data set are in fact unexpected, the 
exposure S in Eqs. (12) and (13) should thus be equal to Sgg. For a fixed number of total 
"observed" events, the larger the background ratio, or, equivalently, the smaller the number of 
real WIMP-induced events, the smaller the required exposure £ — Sgg for accumulating the total 
observed events, and, therefore, the larger the estimated SI WIMP-nucleon coupling. In other 

^°Note that, rather than the mean values, the (bounds on the) reconstructed |/pp are always the median values 

of the simulated results. 

^^As in Refs. [7, 8], in order to avoid complicated calculations of the correlations between the uncertainty on 
estimated by the algorithmic procedure and those on r(Qmin) and Iq, we assumed here that the two data 
sets using the Ge target are independent of each other. 

^^Since {dR/dQ)expt is the measured recoil spectrum before normalized by the exposure £ = Egg, 
-Rexpt(Qmin, Qmax) here is in fact the total number of observed events A/tot = -^sg + ^bg- 




Figure 3: The reconstructed squared SI WIMP-nucleon couplings and the lower and upper 
bounds of the la statistical uncertainties by using mixed data sets from WIMP-induced and 
background events as functions of the input WIMP mass. A ^®Ge nucleus has been chosen as 
the target. The WIMP mass needed in Eq. (13) has been assumed to be known precisely. 
The background ratios shown here are no background (dashed green), 10% (long-dotted blue), 
20% (solid red), and 40% (dash-dotted cyan) background events in the analyzed data sets 
in the experimental energy ranges between and 100 keV. The double-dotted black curve 
indicates the theoretical value of |/pp corresponding to the fixed SI WIMP-nucleon cross section 
(j^p = 10~^ pb. Each data set contains 50 total events on average. Other parameters are as in 
Figs. 1. See the text for further details. 



words, due to extra unexpected background events in our data set, one will use a larger number 
of total events to estimate the SI WIMP coupling, and thus overestimate it. 

More exactly (and mathematically), we can separate the prefactor in the second bracket on 
the right-hand side of Eq. (13) into two terms: 
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It can thus be seen clearly that the prefactor in the expression (13) for estimating |/pp would 
always be overestimated with non-negligible background events. Remind that r((5min) and Iq 
given in Eqs. (10) and (11) are estimated from the measured recoil spectrum {dR/dQ)expt before 
normalized by S (or, equivalently, Ssg). Hence, while the first term on the right-hand side of 
Eq. (23) remains unchanged by increasing the background ratio (and in turn with a decreased 
number of WIMP-induced events), the second term above contributed from residue background 
events causes the overestimate of the reconstructed SI WIMP coupling. Remind also that the 
experimental minimal cut-off energy, Qmin, has been set to be negligible. Thus the first term 
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Figure 4: The products of [rbg/(l — Tbg)] Rsg and (m^ + m^) as functions of for a ^^Ge 
target with different background ratios rbg. The background ratios shown here are 5% (dotted 
magenta), 10% (long-dotted blue), 20% (solid red), and 40% (dash-dotted cyan). Parameters 
are as in Fig. 3. See the text for further details. 



1 /2 

involving Qj^i^ri^g^Qram) in the bracket of the second term above does not contribute to the 
reconstructed |/pp in our simulations shown here; otherwise, the reconstructed |/pp could be 
more strongly overestimated, especially for WIMP masses > 50 GeV (see Figs. 1). 

Moreover, for three cases with background ratios ^ 20% shown in Fig. 3, the larger the input 
WIMP mass, the more strongly overestimated the SI WIMP coupling. However, interestingly, 
once the background ratio rises to > 20% (the dash-dotted cyan curves indicate a background 
ration of 40%), a hump at an input WIMP mass of ~ 60 GeV appears^^. The reason is as 
follows. In the appendix I will show that the second term on the right-hand side of Eq. (23) is 
proportional to the 'WIMP scattering" spectrum {not the "background" spectrum!): 



{Q mm / , T 
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Here {dR/dQ)sg and Rsg are the normalized differential and total event rate of WIMP signals, 
respectively; rbg is the ratio of residue background events in the whole data set. It can be 
understood from Eq. (1) that {dR/dQ)sg and therefore Rsg are functions of the input (true) 
WIMP mass, through not only and rri^^N in the denominator of A defined in Eq. (2), but 
also the transformation constant a in Eq. (5). 

Fig. 4 shows the products of [rbg/(l — rbg)] -Rsg and {m^ + m^) as functions of rri^ for a 
^^Ge target with different background ratios rbg. It can be found that the extra contribution 



-•^■^Remind that the actual values of the "critical" background ratio and the "critical" WIMP mass (with the 
largest systematic deviation) depend in practice strongly on the WIMP scattering spectrum as well as on the 
residue background spectrum and therefore differ from experiment to experiment. 
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Figure 5: As in Fig. 3, except that the WIMP mass in Eq. (13) has been reconstructed by 
using two other data sets with a ^'^Si target and a second '^^Ge target. 

from residue background events, which is proportional to the product of [?"bg/(l — ''"bg)] Rsg and 
{rriy. + mN), has indeed a maximum at a WIMP mass of ~ 75 GeV. Considering the slightly 
decreased |/pp value by increasing the input WIMP mass without background events (the dashed 
green curves in Fig. 3), the total recorded events (including WIMP-induced and background 
events) should thus result in a hump of the reconstructed |/pp at an input WIMP mass of 
~ 60 GeV, once the background fraction in our data set is large enough. 

In Fig. 5 the WIMP mass needed in expression (13) has been reconstructed by using 
two other data sets with a ^®Si target and a second ^^Ge target. As shown in Figs. 2 and 
discussed in Sec. 3.2, due to the contribution from residue background events, if the input 
WIMP mass is light (heavy), the reconstructed mass would be overestimated (underestimated). 
Hence, for input masses < ( ^ ) 150 GeV, the SI WIMP-nucleon couplings reconstructed by 
using three independent data sets would be larger (smaller) than those reconstructed by using 
only one data set with an extra information about the WIMP mass (cf. Fig. 3). In addition, the 
statistical uncertainties on the reconstructed SI WIMP couplings would also be (much) larger. 
However, both Figs. 3 and 5 indicate that one could in principle estimate the SI WIMP-nucleon 
coupling up to a WIMP mass of ~ 1 TeV by using one or three independent data sets with 
maximal 20% background events (solid red). For a WIMP mass of 100 GeV and 20% residue 
background events, the systematic deviation of the reconstructed SI WIMP coupling |/p| (with 
a reconstructed WIMP mass) would in principle be ~ +13% with a statistical uncertainty of 
~lu% ( 3.3%+}^5 for background-free data sets). 

Furthermore, in order to check effects of different background discrimination ability in dif- 
ferent energy ranges and the need of a prior knowledge about an (exact) form of the residue 
background spectrum, in Figs. 6 and 7 we consider two different background windows: between 
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Figure 6: As in Fig. 3, except that the background spectra and/or windows are different: expo- 
nential spectrum between and 50 keV (top), exponential spectrum between 50 and 100 keV 
(middle), and constant spectrum between and 100 keV (bottom). 
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Figure 7: As in Figs. 6, except that the WIMP mass in Eq. (13) has been reconstructed by 
using two other data sets with a ^^Si target and a second '^^Ge target. The background ratios 
shown in two lower frames are no background (dashed green), 5% (dotted magenta), and 10% 
(long-dotted blue). 



and 50 keV and between 50 and 100 keV for the exponential background spectrum as well as 
the rather extrem constant spectrum given in Eq. (20) with a window between and 100 kcV. 

Firstly, from Figs. 6, it can be seen clearly that, by using one data set with up to 20% 
background events and a precisely known WIMP mass as an input information, the SI WIMP- 
nucleon coupling |/p| can in principle be estimated with a maximal ~ +11% systematic deviation 
(for an input WIMP mass of 100 GeV) from the theoretical value and a statistical uncertainty of 
~ ±9%. More importantly, all three cases show almost the same result. This indicates that, once 
the WIMP mass can be known (pretty) precisely, not the exact form of the residue background 
spectrum, but the amount in the analyzed data set could affect (significantly) the reconstructed 
SI WIMP-nucleon couphng. 

In contrast, results shown in Figs. 7 depend strongly on the reconstruction of the WIMP 
mass^^. As discussed in Rcf. [23] and Sec. 3.2, for cases with the exponential background 
spectrum and background windows in the whole experimental possible and low energy ranges, 
the reconstructed WIMP mass could be slightly overestimated (underestimated), once incident 
WIMPs are light (heavy). However, for the case with the exponential spectrum and windows in 
high energy ranges, or the case with the constant spectrum and windows of whole experimental 
possible energy ranges, the reconstructed WIMP mass could be (strongly) overestimated for all 
input WIMP masses. Here the effect of a (strongly) overestimated WIMP mass can be seen 
clearly here. Since for the case shown in Fig. 5 our background spectra are exponential, only 
very few background events could be observed in the energy range between 50 and 100 keV. 
Hence, for the case with background windows only in the low energy ranges (top in Figs. 7), not 
surprisingly, the result of the reconstructed SI WIMP coupling is almost the same as shown in 
Fig. 5. 

However, for the case with the exponential background spectrum and background windows 
in high energy ranges (middle in Figs. 7) or the case with the constant spectrum in whole 
experimental possible energy ranges (bottom), the results are almost the same: the larger the 
background ratio, the more strongly overestimated the SI WIMP coupling, in particular for heavy 
input WIMP masses. Nevertheless, by using (two or) three data sets with background ratios of 
< 10%, one could in principle reconstruct the SI WIMP-nucleon coupling (as well as the WIMP 
mass [23]) pretty well, without knowing the (exact) form of the background spectrum. 

In Figs. 8 we consider a rather light target nucleus: ^^Si. The WIMP mass has been 
assumed to be known precisely (upper) or reconstructed from other two data sets (lower). Only 
the exponential form for residue background spectrum has been considered here. Remind that, 
as found in Refs. [7, 8], with a light target nucleus, e.g.. Si or Ar, the statistical uncertainty on 
the reconstructed SI WIMP-nucleon coupling is larger than that with a heavy nucleus, e.g., Ge 
or Xe. Consequently, for both cases (with a precisely known or a reconstructed WIMP mass), 
the reconstructed SI WIMP couplings as well as the la statistical uncertainties shown in Figs. 8 
are larger than those with '^^Ge shown in Figs. 3 and 5. 

On the other hand, the systematic deviations of the (under) estimated SI WIMP coupling for 
heavy input WIMP masses (m^ > 300 GeV) are smaller for light nuclei than those for heavy 
ones [7, 8]. In addition, as shown in Fig. 9, the background contribution (the second term on 
the right-hand side of Eq. (23)) increases pretty quickly with an increased WIMP mass. Hence, 
for input WIMP masses > 100 GeV, the reconstructed SI WIMP-nucleon couplings with a 
^*Si target nucleus would be more strongly overestimated than those reconstructed with a ^®Ge 
target. Nevertheless, for both cases (with a precisely known or a reconstructed WIMP mass), 

^^Note that in our simulations shown here it was assumed that the spectra and windows of residue background 
events are the same for all three data sets. For practical use with different forms and windows of background 
events in different experiments, one can in principle follow the observations discussed in Sec. 3 and here. 
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Figure 8: As in Figs. 3 (upper) and 5 (lower), except that ^^Si has been used as the (first) target 
for reconstructing /pp. Remind that two data sets with the ^^Si target are different. 
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Figure 9: As in Fig. 4, except that a ^^Si target has been used. 

with ~ 10% - 20% background events in our data sets, the la statistical uncertainty bands could 
in principle still cover the theoretical value of |/p|. For an input WIMP mass of 100 GeV, by 
using three data sets with 10% background events, the systematic deviation of and the statistical 
uncertainty on the reconstructed SI WIMP-nucleon coupling |/p| is ~ +8%~^Iy% (~ "3%^^^^° 
for background-free data sets). 

Finally, considering the progress of detection and background discrimination techniques of 
the next-generation ton-scale detectors, in Figs. 10 we rise the expected number of total events 
in each data set by a factor of 10, i.e., 500 events on average, for both ^^Ge (upper) and ^^Si 
(lower) targets. Here I show only the results with the reconstructed WIMP mass. Since the 
statistical uncertainties shrink now by a factor of ^ 3, the maximal acceptable background ratio 
becomes ~ 5% (i.e., ~ 25 residue background events in each data set). For an input WIMP 
mass of 100 GeV, the systematic deviation would then be ~ +2% (with Ge) and ~ +5% (with 
Si) with a statistical uncertainty of ~ 5% (with Ge) and ~ 7% (with Si). 

5 Summary and conclusions 

In this paper I reexamine the model-independent data analysis method introduced in Refs. [7, 
8] for the estimation of the spin-independent scalar coupling of Weakly Interacting Massive 
Particles on nucleons from data (measured recoil energies) of direct Dark Matter detection 
experiments directly by taking into account a fraction of residue background events, which pass 
all discrimination criteria and then mix with other real WIMP-induced events in the analyzed 
data sets. This method requires neither prior knowledge about the WIMP scattering spectrum 
nor about different possible background spectra; the needed information is the recoil energies 
recorded in direct detection experiments, an unique assumption about the local WIMP density, 
and (occasionally) the mass of incident WIMPs. 

For the mass of incident WIMPs required in this data analysis, we considered two cases: 
known precisely with an overall uncertainty (of 5% of the input WIMP mass in our simulations) 
from other (e.g., collider) experiments as well as reconstructed by using other direct detection 
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Figure 10: As in Fig. 5 and the lower frame of Figs. 8, except that the expected number of total 
events in each data set has been risen by a factor of 10, i.e., 500 events on average, for both the 
™Ge (upper) and the ^^Si (lower) target. The background ratios shown here are no background 
(dashed green), 5% (dotted magenta), and 10% (long-dotted blue) background events in the 
analyzed data sets in the experimental energy ranges between and 100 keV. 



experiments. Our simulations show that, assuming an exponential form for the residue back- 
ground spectrum, with ~ 50 total events in each data set, for both cases the maximal acceptable 
background ratio is ~ 20% (i.e., ~ 10 background events). For a WIMP mass of 100 GeV and 
20% residue background events, the systematic deviation of the reconstructed SI WIMP cou- 
pling |/p| (with a reconstructed WIMP mass) would in principle be ~ -1-13% with a statistical 
uncertainty of ~li4% (~ — 3.3%1}3^° for background-free data sets). 

Furthermore, in order to check effects of different background discrimination ability in dif- 
ferent energy ranges and the need of a prior knowledge about an (exact) form of the residue 
background spectrum, we considered two different background windows: between and 50 keV 
and between 50 and 100 keV for the exponential background spectrum as well as the rather 
extrem constant spectrum with a window between and 100 keV. It has been found that, with 
a precisely known WIMP mass, all three cases show almost the same result as that for the ex- 
ponential background spectrum with a window of the whole experimental possible energy range. 
This indicates that, once the WIMP mass can be known (pretty) precisely, not the exact form 
of the residue background spectrum, but their amount in the analyzed data set could affect 
(significantly) the reconstructed SI WIMP-nucleon coupling. 

On the other hand, we considered a rather light target nucleus: ^^Si with both a precisely 
known and a reconstructed WIMP mass. For both cases the reconstructed SI WIMP couplings 
would be more strongly overestimated and the la statistical uncertainties would also be larger 
than those with ^^Ge. Nevertheless, with ~ 10% - 20% background events in our data sets, the 
1(7 statistical uncertainty bands could in principle still cover the theoretical value of |/p|. For 
an input WIMP mass of 100 GeV. by using three data sets with 10% background events, the 
systematic deviation of and the statistical uncertainty on the reconstructed SI WIMP-nucleon 
coupling |/p| is +8%-i7% ("^ ~3%+^7% for background-free data sets). 

Finally, for rather next-generation ton-scale detectors, we considered the use of data sets of 
C(500) events for both ^^Ge and ^^Si targets. Our results show that, with a maximal background 
ratio of 5% (i.e., ~ 25 total events in each data set), one could in principle still reconstructed the 
SI WIMP-nucleon coupling pretty well: for an input WIMP mass of 100 GeV, the systematic 
deviation would be ~ +2% (with Ge) and ~ -|-5% (with Si) with a statistical uncertainty of 
~ 5% (with Ge) and - 7% (with Si). 

In summary, as the third part of the study of the effects of residue background events in di- 
rect Dark Matter detection experiments, we considered the estimation of the SI WIMP-nucleon 
coupling. Our results show that, with currently running and projected experiments using detec- 
tors with 10^^ to 10^^^ pb sensitivities [29, 20, 30, 31] and < 10^'' background rejection ability 
[19, 21, 22, 18], once one or more experiments with different target nuclei could accumulate a 
few tens events (in one experiment), we could in principle already estimate the SI coupling of 
Dark Matter particles on ordinary matter with a reasonable precession, or at least give an upper 
bound on that, even though there could be some background events mixed in our data sets for 
the analysis and the reconstructed value would thus be overestimated. Moreover, although two 
forms for background spectrum and three windows for residue background events considered in 
this work are rather naive; one should be able to extend our observations/discussions to predict 
the effects of possible background events in their own experiment. Hopefully, this will encourage 
our experimental colleagues to present their (future) results as the "most possible area(s)" in the 
parameter space of different extensions of the Standard Model of particle physics and in turn to 
offer stringenter information for identifying (WIMP) Dark Matter particles at colliders as well 
as for predicting spectra in indirect Dark Matter detection experiments. 
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A Formulae needed in Sec. 2 

Here 1 list all formulae needed for the model-independent data analysis procedure used in Sec. 2. 
Detailed derivations and discussions can be found in Refs. [12, 6]. 



A.l Estimating r(Qmin) and /n(Qmin, Qmax) 

First, consider experimental data described by 
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Here the total energy range between Qmin and Qmax has been divided into B bins with central 
points Qn and widths 6„. In each bin, Nn events will be recorded. Since the recoil spectrum 

dR/dQ is expected to be approximately exponential, the following ansatz for the measured ieco][ 
spectrum [before normalized by the experimental exposure £) in the nth bin has been introduced 
[12]: 
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Here r„ is the standard estimator for (dR/ dQ)e^pt Q — Qn- 



bn ' 



(A3) 



kn is the logarithmic slope of the recoil spectrum in the nth Q— bin, which can be computed 
numerically from the average value of the measured recoil energies in this bin: 
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The error on the logarithmic slope kn can be estimated from Eq. (A4) directly as 
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Qs,n in the ansatz (A2) is the shifted point at which the leading systematic error due to the 
ansatz is minimal [12], 
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Note that Qs,n differs from the central point of the nth bin, Q^. Prom the ansatz (A2), the 
counting rate at Q — Qmm can be calculated by 



and its statistical error can be expressed as 
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Finally, since all are determined from the same data, they are correlated with 
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where the sum runs over all events with recoil energy between Qmin and Qmax- And the correlation 
between the errors on r{Q-[am), which is calculated entirely from the events in the first bin, and 
on /„ is given by 
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note that the sums here only count in the first bin, which ends at Q = Qmin + ^i- 

On the other hand, with a functional form of the recoil spectrum (e.g., fitted to experimental 
data), {dR/ dQ)e:x.pt, one can use the following integral forms to replace the summations given 
above. Firstly, the average Q— value in the nth bin defined in Eq. (A5) can be calculated by 
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For /n(Qmm, Qmax) givcn in Eq. (11), we have 
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and similarly for the covariance matrix for 7„ in Eq. (A12), 
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Remind that (di?/d(5)expt is the measured recoil spectrum before normahzed by the exposure. 
Finally, /^(Qmin, Qmin + ^'i) needed in Eq. (A13) can be calculated by 
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Note that, firstly, r(Qmm) and IniQmim Qmm + bi) should be estimated by Eqs. (A9) and (A17) 
with ri, ki and estimated by Eqs. (A3), (A4), and (A8) in order to use the other formulae for 
estimating the (correlations between the) statistical errors without any modification. Secondly, 
'"(Qmin) and /n((5min, Qmax) estimated from a scattering spectrum fitted to experimental data 
are usually not model-independent any more. 



A. 2 Determining the WIMP mass m 



X 



By requiring that the values of a given moment of fi{v) estimated by Eq. (9) from two exper- 
iments with different target nuclei, X and Y, agree, appearing in the prefactor a" on the 
right-hand side of Eq. (9) can be solved analytically as [5, 6]: 



m 



with TZn,{x,Y) defined by 



(A18) 
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,x) + (n + l)/„,x 



2Qrnin X^^(Qmm,x)/-^l(Qmm,x) + Iq,X 



1/n 



(A19) 



and TZnY can be defined analogously. Here n ^ m(x.Y) and F(^x,Y)iQ) masses and the 

form factors of the nucleus X and Y, respectively, and r(x,y)(Qmm,(x,y)) refer to the counting 
rates for the target X and Y at the respective lowest recoil energies included in the analysis. Note 
that the general expression (A18) can be used either for spin-independent or for spin-dependent 
scattering, one only needs to choose different form factors under different assumptions; the form 
factors needed for estimating In,{x,Y) by Eq- (H) or (A15) are thus also different. 

By using the standard Gaussian error propagation, a lengthy expression for the statistical 
uncertainty on m^\(^^„-^ given in Eq. (A18) can be obtained as 
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Here a short-hand notation for the six quantities on which the estimate of depends has been 
introduced: 



Cl,X 



in,X 



C2,X = -^0,. 



C3,X = rx{Qmin,x) ] 



(A21) 



and similarly for the Q^y. Estimators for cov(cj,Cj) have been given in Eqs. (A12) and (A13). 
Explicit expressions for the derivatives of TZn^x with respect to q^x are: 
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(A22a) 
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and 
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(A22c) 



explicit expressions for the derivatives of 7^„,y with respect to Q^y can be given analogously. Note 
that, firstly, factors 7^n,(x,y) appear in all these expressions, which can practically be cancelled 
by the prefactors in the bracket in Eq. (A20). Secondly, all the Iq,{x,y) In,{x,Y) should be 
understood to be computed according to Eq. (11) or (A15) with integration limits Qmin and 
Qmax specific for that target. 

On the other hand, since |/pp in Eq. (13) is identical for different targets, it leads to a second 
expression for determining [6]: 



m 
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Here m(^x,Y) ^(x,y) has been assumed, and TZa,{x,Y) defined by 
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and similarly for TZcr,Y'i ^{x,y) here arc the experimental exposures with the target X and Y. 
Similar to the analogy between Eqs. (A18) and (A23), the statistical uncertainty on m^|^ given 
in Eq. (A23) can be expressed as 
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where I have used again the short-hand notation in Eq. (A21); note that ci (x,y) = In,{x,Y) '^^ 
not appear here. Expressions for the derivatives of TZ^^x can be computed from Eq. (A24) as 
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(A26b) 



and similarly for the derivatives of 7^o-,y- Remind that factors ncr,{x,Y) appearing here can also 
be cancelled by the prefactors in the bracket in Eq. (A25). 



In order to yield the best-fit WIMP mass as well as to minimize its statistical uncertainty 
by combining the estimators for different n in Eq. (A18) with each other and with the estimator 
in Eq. (A23), a function has been introduced as [6] 
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the other nmax + 2 functions fi^y can be defined analogously. Here rimax determines the highest 
moment of fi{v) that is included in the fit. The are normalized such that they are dimen- 
sionless and very roughly of order unity in order to alleviate numerical problems associated with 
the inversion of their covariance matrix. Note that the first Umax + 1 fit functions depend on 
only through the overall factor a and in Eqs. (A28a) and (A28b) is now a fit parameter, 
which may differ from the true value of the WIMP mass. Finally, C in Eq. (A27) is the total 
covariance matrix. Since the X and Y quantities are statistically completely independent, C can 
be written as a sum of two terms: 



dj = COV fj,x) + GOV {fi^Y, fjx) 



(A29) 



The entries of the C matrix given here involving basically only the moments of the WIMP 
velocity distribution can be read off Eq. (82) of Ref. [12], with an slight modification due to the 
normalization factor in Eq. (A28a)^^: 



GOV (/i, /,) = 



fi fj cov(/o, /o) + + + l)cov(/,, /,) 

- a^{j + l)fi cov(7o, Ij) - oc\i + l)fj cov(7o, h) 

+ A^jCT2(r(Q„,in)) - (A/j + Djfi) cov(r(Qmin), h) 

+ &{3 + 1) A cov(r(Qmin), Ij) + oi\i + l)Dj cov(r(Qmin), h) 



Here I used 
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^^Since the last fi defined in Eq. (A28b) can be computed from the same basic quantities, i.e., the counting 
and the integrals Iq, it can directly be included in the covariance matrix. 
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Finally, since the basic requirement of the expressions for determining given in Eqs. (A18) 
and (A23) is that, from two experiments with different target nuclei, the values of a given moment 
of the WIMP velocity distribution estimated by Eq. (9) should agree, the upper cuts on fi{v) 
in two data sets should be (approximately) equal^^. Since fcut = a^Qmax, it requires that [6] 
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Note that a defined in Eq. (5) is a function of the true WIMP mass. Thus this relation for 
matching optimal cut-off energies can be used only if m^^ is already known. One possibility to 
overcome this problem is to fix the cut-off energy of the experiment with the heavier target, 
minimize the x^(m^) function defined in Eq. (A27), and then estimate the cut-off energy for 
the lighter nucleus by Eq. (A34) algorithmically [6]. 



A. 3 Statistical uncertainty on |/pp 

By using the standard Gaussian error propagation, the statistical uncertainty on 
by Eq. (13) can be given as 
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Here the statistical error on 1/Nm. can be given from Eq. (A31) directly as 
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For the case that one has only two data sets with different target nuclei, X and F, one of these 
two data sets will then be needed for reconstructing the WIMP mass rriy^ and also for estimating 
in Eq. (13). The uncertainties on and 1/jVm are thus correlated. Assuming that the 
WIMP mass is reconstructed by Eq. (A18), and target X{Y) is used for estimating 1/jVm, the 
covariance of ?7t.^|^^„^ and l/A/'m,(x,y) can be obtained by modifying Eq. (A20) slightly as 
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Here the threshold energies have been assumed to be negUgible. 
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For the case that the WIMP mass is reconstructed by Eq. (A23), one can also modify Eq. (A25) 
to obtain that 
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Note that, firstly, in the above expressions we have to use (mx — my) instead of |mx — my| in 
Eqs. (A20) and (A25); for expressions with the Y target, there is an additional "— (minus)" 
sign. Secondly, the algorithmic process for matching the experimental maximal cut-off energies 
of two experiments used for the reconstruction of the WIMP mass can also be used with the 
basic expressions (A18) and (A23). For this case and the lighter nucleus is used for estimating 
l/-^m, the energy range of the sum in Eq. (A12) or of the integral in Eq. (A16) as the estimator 
for the covariance of /„ should be modified to be between Qmin and the reduced maximal cut-off 
energy of the lighter nucleus. 



B Proportionality of r'bg(Qmin) and /o,bg to R^^ 

The spectrum of residue background events before normalized by the experimental exposure 
8 — Ssg can be expressed as 



'dR\ 



a£.,[^] . (A39) 

bg,expt 

Here a is a proportional constant, £sg is the required exposure to observe the expected "WIMP 
signal" (not total) events, which can be estimated theoretically by 

^^^j dQ , (A40) 
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where A^tot and A'gg are the number of the total and WIMP signal events, respectively; < rbg < 1 
is the ratio of residue background events in the whole data set. On the other hand, the number 
of residue background events in the data set can be given by 
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{dR/dQ)bg in Eq. (A39) and here is a (simphfied) analytic form of the background spectrum, 
e.g., {dR/dQ)^g^e^ in Eq. (21) and ((ii?/o?(5)bg,const in Eq. (20). 

Similar to Eq. (A15), /„bg can be estimated from the background spectrum (dit!/d(5)bgexpt 

by 
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Hence, the second term involving rbg((5min) and /o,bg(Qmin, Qmax) on the right-hand side of 
Eq. (23) can be given as 
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Finally, by combining Eqs. (A40) and (A41), the proportional constant a can be calculated by 
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Remind that, while the signal spectrum {dR/dQ)sg is a function of the WIMP mass m^, the 
background spectrum {dR/dQ)\yg should in general be independent oi m^. 
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